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The gravitational wave (GW) signals emitted by a network of cosmic strings are reexamined in 
view of the possible formation of a network of cosmic superstrings at the end of brane inflation. The 
reconnection probability p of intersecting fundamental or Dirichlet strings might be much smaller 
than 1, and the properties of the resulting string network may differ significantly from those of 
ordinary strings (which have p = 1). In addition, it has been recently suggested that the typical 
length of newly formed loops may differ by a factor e <C 1 from its standard estimate. Here, we 

■ analyze the effects of the two parameters p and e on the GW signatures of strings. We consider both 
the GW bursts emitted from cusps of oscillating string loops, which have been suggested as candidate 
sources for the LIGO/VIRGO and LISA interferometers, and the stochastic GW background, which 
may be detectable by pulsar timing observations. In both cases we find that previously obtained 
results are quite robust, at least when the loop sizes are not suppressed by many orders of magnitude 
relative to the standard scenario. We urge pulsar observers to reanalyze a recently obtained 17- year 

P_h ' combined data set to see whether the large scatter exhibited by a fraction of the data might be due 

I , to a transient GW burst activity of some sort, e.g. to a near cusp event. 

(N ! I. INTRODUCTION 

> ; 

Cosmic strings can be formed as linear defects at a symmetry breaking phase transition in the early universe 
and can give rise to a variety of observable phenomena at the present cosmic age. String formation, evolution, and 
observational effects have been extensively studied in the 1980's and 90's (for a review see [1,2]). This exploration 
t-H ' was to a large degree motivated by the string scenario of structure formation [3,4], which requires strings of the grand 
unification energy scale, rj ~ 10 16 GeV, or Gp ~ (rj/Mp) 2 ~ 10~ 6 . ( Here, G is Newton's constant, p is the string 
tension and M p ~ 10 19 GeV is the Planck mass. Gp is a dimensionless parameter characterizing the gravitational 
interactions of strings.) This scenario is now disfavored by the CMB observations, but strings of a somewhat lower 

■ energy scale would be consistent with the data, and their detection would of course be of great interest. The current 
O i| CMB bound on strings is [5,6] Gp < 6.1 x 10 -7 . There is also a case of potential string detection with a similar value 

of Gp. Two nearly identical galaxies are observed at angular separation of 1.9 arc sec, suggesting gravitational lensing 
' by a cosmic string with Gp ~4x 10~ 7 [7]. 

, Cosmic strings can also be detected through the gravitational wave (GW) background produced by oscillating 
string loops [8]. This background, which ranges over many decades in frequency, has been extensively discussed in the 
literature [8-13]. The analysis of eight years of millisecond pulsar timing observations has led to setting rather stringent 
(95% confidence level) limits on the GW contribution to the cosmological closure density Q g = p g /p c = S^Gpg/ (3Hq): 
ilgh 2 < 6 x 1CP 8 according to the original analysis [14], or fl g h 2 < 9.3 x 10~ 8 according to the Bayesian approach of 
[15]. (Here, h = _ff /(100km/s/Mpc). Note that h 2 ~ (65/100) 2 ~ 0.42.) As we shall review below, the corresponding 
bound on Gp is c 3 ' 2 Gp < 10 -7 , where c denotes the (mean) number of cusp events per oscillation period of a string 
loop. 

Until recently, it appeared that the gravitational effects of strings with Gp <C 10~ 7 are too weak to be observable. 
However, it has been shown in [16,17] (following a suggestion in [18]) that GW bursts emitted from cusps of oscillating 
loops should be detectable by LIGO and LISA interferometers for values of Gp as low as 10 -13 . 

During the last few years, there have been several important developments that motivate us to reexamine the 
GW signatures of strings. First, there has been a renewed interest in the possibility [19] that fundamental strings 
of superstring theory may have astronomical sizes and play the role of cosmic strings. In particular, it has been 
argued [20-22] that fundamental (F) and D-string networks can naturally be formed at the end of brane inflation. 
In this scenario [23], inflation is driven by the attractive potential between a D-brane and an anti-D-brane, and 
strings are produced when the branes eventually collide and annihilate. The rather stringent requirements allowing 
for the stability of such cosmic superstrings can be met in some scenarios [22,24,25]. The predicted string tensions 
are [20,21,26,24] 10 -11 < Gp < 10~ 6 and appear to be comfortably within the range of detectability by LIGO and 
LISA. However, the analysis of GW bursts in [16,17] may not be directly applicable in this case, since the properties 
and evolution of F and D-string networks may differ in significant ways from those of "ordinary" cosmic strings. 
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If both F and D-strings are produced, they can form an interconnected network, in which F and D-strings join and 
separate at 3-way junctions [24,22]. Each junction joins an F-string, a D-string, and a (l,l)-string, which is a bound 
state of F and D. (p, ^-strings, which are bound states of p F-strings and q D-strings with p, q > 1 can also be formed. 
The evolution of FD-networks is similar to that of monopole-string Z 3 networks, in which each monopole is attached 
to three strings. Simulations of Z 3 network evolution suggest [27,28] that the typical distance between the monopolcs 
scales with the cosmic time, L ~ -yt, where the coefficient 7 depends on the rate of energy loss by the network. But 
if the main energy loss mechanism is gravitational radiation, as the case may be for an FD network, then it has been 
argued in [27] that the energy dissipation by the network is rather inefficient, so it quickly comes to dominate the 
universe. If this picture is correct, then models predicting FD networks are ruled out. It is conceivable, however, that 
the main energy loss mechanism of FD networks is not GW emission, but chopping off of small nets, similar to closed 
loop production by ordinary strings. The negative verdict on this type of models can then be avoided. This issue can 
only be resolved with the aid of new, high-resolution numerical simulations of FD networks. 

In this paper we shall focus on models where only one type of string is formed. Still, the string evolution may differ 
from that of ordinary strings, because the reconnection probability p for intersecting strings may be significantly 
smaller than 1. When ordinary strings intersect, they always reconnect [29,30]. For intersecting F-strings, the 
reconnection probability is suppressed by the string coupling, gl < 1. Moreover, strings moving in a higher-dimensional 
bulk can avoid intersection much more easily than strings in 3 dimensions [22,21]. The string propagation in the bulk 
is expected to be restricted by bulk potentials, but the effective reconnection probability can still be reduced by an 
order of magnitude or so. Analysis in [25] suggests reconnection probabilities in the range 

10~ 3 <P<1 (1-1) 

for F-strings and 

0.1<P<1 (1-2) 

for D-strings. We thus need to analyze the effect of a reduced reconnection probability on the GW burst statistics 
and on the stochastic GW background. 

Another interesting recent development has been the analysis by Siemens and Olum [31] of the gravitational 
radiation from counter-streaming wiggles on long strings. They showed that this radiation is much less efficient 
in damping the small-scale wiggles than originally thought. This may result in much smaller sizes of closed loops 
produced by the network [32] than previously assumed. This effect can be quantified by the dimensionless parameter 
e defined in Eq. (2.9) below. Refs. [16,17] had assumed (besides p — 1), the "standard" value e = 1, and we need to 
see how a different value of e might affect the GW burst statistics. 

Last, but certainly not least on our list of recent developments is the potential improvement in the sensitivity to a 
GW background of pulsar timing observations. Indeed, these have been recently extended to a 17- year data set [34]. 
We shall discuss below what kind of limits on G[i follow from this extended data set. 

The main goal of the present paper is to analyze how the amplitude and frequency of GW bursts from strings 
and the intensity of the stochastic GW background depend on the parameter e measuring the characteristic size of 
closed loops, and on the string reconnection probability p. The paper is organized as follows. In the next Section we 
outline the relevant features of string evolution. GW bursts from cusps are discussed in Section III. The stochastic 
background is discussed in Section IV, where we also discuss bounds on G[i based on millisecond pulsar observations. 
Our conclusions are summarized in Section V. 



II. STRING EVOLUTION 
A. Standard scenario 

An evolving string network consists of two components: long strings and sub-horizon closed loops. The long 
string component is characterized by the following parameters: the coherence length £(t), defined as the distance 
beyond which the directions along the string are uncorrelated, the average distance between the strings L(t), and the 
characteristic wavelength of the smallest wiggles on long strings, l w iggles(t)- The standard picture of cosmic string 
evolution (based on the assumptions p = 1 and e = 1), asserts that (in units where the velocity of light is set to one) 

L rt (t)~f t (t)~t, (2.1) 

and 
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lwiggles{t) ~ at, (2.2) 

where a is a constant whose value is specified below. 

Strings move at relativistic speeds, and each long string intersects itself or another long string about once every 
Hubble time t. As a result, one or few large loops of size L st (t) ~ t are produced per Hubble time, which then shatter, 
through multiple sclf-intcrscctions, into a large number of small loops, whose size is comparable to the wavelength of 
the wiggles (2.2). In other words, when loops are just formed they have a typical size 

l(t) ~ at. (2.3) 

This equation, referring to the typical size of just formed loops, will hold throughout this paper (including in the "non 
standard" cases discussed below), and defines the meaning of the dimensionless loop-length parameter a. 

Note that in order for the characteristic length of the network to scale with cosmic time as in (2.1), it is necessary 
for long strings to discharge a sizeable fraction of their length (~ t per Hubble volume per Hubble time) in the form 
of closed loops. In other words, the number of loops formed per Hubble volume per Hubble time, say Ni, is on the 
order of 

(2.4) 

The loops oscillate and lose their energy by gravitational radiation at the rate 

d£/dt~TGn 2 , (2.5) 

where r ~ 50 [1] is a numerical coefficient. The lifetime of a loop of length l(t) ~ at and energy £ ~ ~ fiat is 

t - (a/TGfi)t. (2.6) 

The standard scenario assumes that the value of a is determined by the gravitational damping of small-scale wiggles. 
With the naive, old estimate of the damping, one finds [11] 

a st ~ TG/j,. (2.7) 

Then the lifetime of loops formed at time t is t ~ t, so that, at any moment, the number of loops per Hubble volume 
is the same as the number of loops formed per Hubble volume per Hubble time, as given by Eq. (2.4). Therefore the 
number density of loops in the standard scenario is 

n^t)-^ 1 *- 3 -^)- 1 *- 3 . (2.8) 

In the following subsections we shall consider the modifications to the result Eq. (2.8) of the standard scenario brought 
by two possible effects: (i) a small reconnection probability p <C 1, and (ii) a value of the loop- length parameter a 
differing from the standard value Eq. (2.7) by being cither smaller, or larger than it. We can quantify the effect (ii) 
with a new parameter 

e = a/a st = a/TGfj,. (2.9) 

Note again that we shall always define a by Eq. (2.3) (concerning the typical length of newly formed loops), which 
therefore holds true both in the standard scenario, and in the extended scenarios studied here. On the other hand, 
Eq. (2.7) will cease to hold in the extended scenarios. 



B. Small reconnection probability: p<l 

On sub-horizon scales, the strings are straightened out by the expansion of the universe, and thus we expect, as 
before, that the string coherence length is 

£(t)~t. (2.10) 

If the reconnection probability is p <C 1 , then one intersection per Hubble time is not sufficient for ensuring the scaling 
of long strings. In order to have one reconnection, a long string needs to have ~ p~ x intersections per Hubble time. 
This means that the number of such strings per Hubble volume should be ~ p^ 1 . 
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The typical inter-string distance L(t) can be estimated from comparing two different estimates of the mean energy 
density in long strings, say pi on g- On the one hand, considering that a Hubble volume ~ t 3 contains ~ p~ x long 
strings yields pi on g ~ p^ x [itjt 3 = p/(pt 2 ). On the other hand, the definition of L is that there should be ~ 1 string 
segment of length L in a spatial volume ~ L 3 , so that pi ong ~ fiL/ 'L 3 = p/L 2 . Equating the two estimates yields 

L(t)~p 1/2 t, (2.11) 

instead of Eq. (2.1). 

There have been some conflicting claims in the recent literature regarding the values of L(t) and in the regime 
of p <C 1. Jones, Stoica and Tye [21] find that L(t) ~ £(t) ~ pf. We note, however, that they derive this result from 
the one-scale model [35], which assumes from the outset that L(t) ~ £(t). This assumption is not necessarily justified 
when p <C 1. Sakellariadou [36] argues that £(t) ~ £(t) ~ p x / 2 t. Her analysis is based on numerical simulations 
of string evolution in flat spacetime. Eq. (2.11) was indeed originally obtained in such simulations [37], but we are 
not aware of any results indicating that ~ L(t) in this regime. On the contrary, preliminary results of a new, 
high- resolution flat-space simulation indicate that is substantially larger than L(t) [38]. 

A definitive picture of string evolution with low intercommuting probability can only be obtained from numerical 
simulations in an expanding universe. Flat-space simulations do not include the important effect of string stretching 
on sub-horizon scales. Simple analytic estimates do not account for the possible effects of the wiggles. For example, 
colliding wiggly strings may intersect at more than one point, thus increasing the effective reconnection probability. 
If the wiggliness of strings is increased at low p, their velocity will be reduced, leading to a further decrease of L(t). 
String simulations with p< 1 will hopefully be performed in the near future. For the time being, we shall assume 
that Eqs.(2.10), (2.11) give a reasonable approximation and use them in the rest of the paper. 

As before, a significant fraction of the total string length within a Hubble volume ~ p~ x t should go into loops each 
Hubble time, so that the number of loops formed per Hubble volume per Hubble time is now 

lit) pa 

i.e p~ x larger than the standard result (2.4). In order to estimate the corresponding loop density at any moment, one 
must take into account the lifetime of these loops (which depends on the value of e). This will be discussed in the 
following subsections. 



C. Small loops: a <C TGp, i.e. e < 1 

As we already mentioned, recent analysis in [31] has shown that the gravitational radiation from counter-streaming 
wiggles on long strings is much less efficient in damping the wiggles than originally thought. If indeed a is determined 
by gravitational back-reaction, then the new analysis shows [32] that its value is sensitive to the spectrum of small- 
scale wiggles, and is generally much smaller than TGp. A simple model for the spectrum of wiggles introduced in [32] 
yields 

a ~ (TGfj,) n (2.13) 

with n = 3/2 in the radiation era and n = 5/2 in the matter era. In other words, the quantity (2.9) would be of 
order e ~ (TGp) m , with m = n — 1, i.e. m = 1/2 in the radiation era and m — 3/2 in the matter era. As TG/i 
is observationally restricted to be TG/j, < 1CP 5 , and might turn out to be -C 10~ 5 , this gives very small values of 
a < 3 x 10~ 8 and a<3x 10~ 13 , and e<3x 1CP 3 and e<3x 1CP 8 , for radiation and matter eras, respectively. In 
view of such possible drastic changes in orders of magnitude, it is important to assess the effect of e <C 1 on string 
network statistics and the corresponding GW observables. 

As said above, for a -C TG/i, the number of loops produced per Hubble time per Hubble volume is given by 
Eq. (2.12). However, to compute the number of loops at any moment we must now take into account the fact that 
the lifetime of the loops (2.6) is r <C t. Therefore, only a small fraction of these loops will be present at any given 
time, N' ~ (j/t)N ~ l/(pTGp). The corresponding loop density is therefore 

Note that, compared to the standard result Eq. (2.8), the value of n(t) in scenarios extended by the two parameters 
p <C 1 and e <C 1 is independent of e, and exhibits a simple dependence cx upon p. 
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D. Large loops: a ^> TG/i, i.e. e S> 1 



As mentionned above, the possibility a <C TG^i, i.e. e ^ 1 will be our main focus in this paper, because of the 
recent findings of [31,32]. However, the opposite case a 3> FGfi, i.e. e>l cannot, at this stage, be dismissed. Indeed, 
the spectrum of the wiggles is expected to be a power law decaying towards shorter wavelengths [32] . This raises the 
possibility that for sufficiently small wavelengths the wiggles may be too small to have an effect on loop formation. 
The size of the loops will then be determined by the dynamics of the network, and gravitational back-reaction will 
play no role. 

The possibility of a ^> TG/j has been discussed in [1]. We recall that the parameter a is defined by requiring that 
the typical length of newly formed loops is ~ atf, where tf denotes the time of formation. When a S> TGfi, such 
loops will survive over many Hubble times, and the loops extant at any given cosmological time t will be obtained 
by integrating over the loops formed on all formation times tf < t. One finds that there is a distribution of loops 
with sizes in the range < I < at, with the dominant contribution to the number density, and to the GW burst rate, 
coming, at cosmological time t, from loops of typical size 

I ~ TG^it. (2.15) 

These dominant loops at time t were formed at the parametrically smaller typical time tf ~ (TGfj,/a)t <C t. Their 
density at that (formation) time was 



~ (2 - 16) 



and the density at time t is 

n(t) 



<tf) 



a(t) 
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n f , (2.17) 



where a(t) denotes the cosmological scale factor. If we consider loops formed in the matter era, tf > t eq (which will 
indeed be the most important for GW observations), the factor (a(t /) / a(t)) 3 — {tf/t) 2 ~ (TGfi/a) 2 , so that the loop 
number density is 

n <*> ~ wh- (2 - 18) 

It is interesting to note that the final result (2.18) for the loop density (in the matter era) when e ^> 1 coincides with 
the result (2.14) obtained in the opposite case e <C 1, and that both results are independent of e. Note, however, that 
the typical size of the loops at time t are different. In the case where e <C 1 this typical size is l(t) ~ at = eTGfit, 
while in the case e » 1 it is l(t) ~ TG^it. In other words, the case e > 1 can be effectively treated by taking the limit 
e — > 1 of the case e < 1 (while keeping the effect of p) . Another way to say this is to introduce the notion of effective 
loop- length parameter a e f / , defined by writing that the typical size of the loops which dominate the loop density at 
cosmic time t is 

l(t) ~ aef ft. (2.19) 

Note that Eq. (2.19) refers to the typical size of loops surviving at some cosmic time t, while Eq. (2.3) referred to the 
typical size of newly formed loops. Correspondingly, we can define e e // = a e ff/TGfj,. With this definition, one has 
a e ff = a when a < YG[i (i.e. e e ff = e when e < 1), and a e ff — TGfi when a > TGfi (i.e. e e ff = 1 when e > 1). 
Note that e e ff is never greater than 1. [One could approximately write the link e e ff = e/(l + e).] 



E. This paper 

The bottom line is that the value of a is presently unknown. Numerical simulations of string evolution [11,39] give 
loops that are too small to be resolved, so only an upper bound on a can be obtained, 

a < 1(T 3 . (2.20) 

The next generation of string simulations, which are now being developed, are expected to improve this bound 
considerably In this paper, we shall use a, or equivalently e defined by Eq. (2.9), as a free parameter and will allow 
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it to take values smaller, as well as greater than its standard value a s t = TG/i, corresponding to e st = 1. Similarly, we 
shall treat the reconnection probability p as a free parameter and consider cases of both high (p = 1) and low (pC 1) 
reconnection probability. 

It should be emphasized though that the parameters p and a do not appear in the theory on equal footing. The 
reconnection probability p and the string tension /i arc true parameters, in the sense that they can take different 
values in different cosmic string models (e.g., fundamental strings, D-strings, or "ordinary" strings). On the other 
hand, the parameter a only reflects incompleteness of our understanding of string evolution and will eventually be 
determined, possibly as a function of G[i and p. 1 

III. GRAVITATIONAL WAVE BURSTS 

The main result of the GW burst analysis in Refs. [16,17] is the expression for the typical dimensionless amplitude 
of cusp-generated bursts, observed in an octave of frequency around frequency /, that one can expect to detect at a 
given occurrence rate N (say, one per year), 

M/) ~ G M a 2 /3 (/fo) -i/3 fl[y (oid) ]e(1 _ dm[aj L Zm(2 /oid) )D . (31) 
Here, to ^ 2/(3H ) ~ 1.0 x 10 10 yr ~ 10 7 5 s is the present cosmic time, 

y<?*)(NJ) = 10- 2 (iV/ C )i a 8/3 (/to) 2/3 , (3.2) 

c is the average number of cusps per period of loop oscillation, and the superscript (old) refers to the fact that Eq.(3.2) 
was derived using the old (standard) string evolution scenario with a ~ TG^i. The function g[y] in (3.1) is given by 

g[y] = y-V3(i + j,)-W33 (1 + y /y eq f/i\ (3 . 3) 

This is an interpolating function which represents the power-law behavior of in three different regimes: y < 1, 

1 < y < y eq , and y > y eqi where y eq — z\ q and z eq ~ 10 3 ' 94 is the rcdshift of equal matter and radiation densities. 
The three regimes correspond to loops radiating, respectively, at z < 1, 1 < z < z eq , and z > z eq . Finally, the last 
factor 9 in Eq.(3.9) involves Heaviside's step function 0, with the argument 1 — 8 m obtained by inserting Eq.(3.2) 
into the function 

z m (y) = y"\l + y yi*\l + yly eq )-^\ (3.4) 
and then by inserting the result into the function 

m (a, /. z ) = (afto)- 1/3 (l + ^) 1/6 (1 + z/ze q ) 1,& . (3.5) 

Physically, the quantity 9 m (/) is related to the (integer) mode number m of the Fourier decomposition of the grav- 
itational radiation emitted by a loop by \m\ ~ (0 m (f))~ 3 - Note also the link \m\ <~ (1 + z)fl [17], where I denotes 
the length of the loop. The 6-function factor serves the purpose of restricting the burst signals to the values m < 1, 
corresponding to |m| > 1. In view of Eq. (3.5), this step function is equal to one when the product afto is larger 
than 1 and z < 1. This case covers many cases of physical interest, especially when considering GW frequencies in 
the LIGO or LISA bands. This is why, in most of our analytical discussion below we shall only consider the other 
factors in Eq. (3.1). However, we shall include the effect of this cut-off factor in our plots below, and we shall see 
that is plays a crucial role for the GW frequencies of relevance in pulsar timing observations, and that it also starts 
playing an important role for LIGO and LISA signals, when e gets much smaller than 1. 

The log-log plot of as a function of G/i, calculated for the standard scenario, and in absence of the cut-off 
brought by the ©-function factor, is made up of three straight lines representing the three different regimes mentioned 
above. Now we would like to find out how these lines are modified when e ^ 1 and/or p < 1. By looking at the 
derivation of Eqs.(3.1)-(3.3) in Refs. [16,17] one finds that the explicit occurrences of the notation a concerned two 
different aspects of a string network: (i) either a was used to parametrize the typical size of a loop, in the sense of 
Eq. (2.3), or better of Eq. (2.19) as one is interested in the typical size of a loop at cosmic time t, or, (ii) a entered as 



1 The parameter c, the average number of cusps per loop, which is introduced in the next Section, is also not a true parameter 
and will hopefully be determined from numerical simulations. 



6 



a factor in the loop density, written as n(t) ~ a _1 i~ 3 . The first usage of a is consistent with the definition Eq. (2.3) 
of a in the extended scenarios considered here in the most interesting case where e < 1; as said above, we shall treat 
the opposite case by replacing a — ► a e ff, i.e. by taking the e — > 1 limit. But then, if a is so defined, the second usage 
must be corrected for because the loop density in extended scenarios is not given by n(t) ~ a -1 ^ 3 , but rather by 
Eqs. (2.14) and (2.17) (the latter becoming Eq. (2.18), which is the same as (2.14) in the matter era case). In other 
words, we just need to correct (at least when e < 1), the loop density used in Rcfs. [16,17] by the factor a/pTGfi = e/p. 
An easy way to accomplish this is to note that the modified loop density can be accounted for by adjusting the value 
of the number of cusp events per loop oscillation, c. Specifically, for e < 1 we need to make a replacement 

c--. (3.6) 
P 

Alternatively, as c always enters the final results in the combination N/c, we could account for the modified loop 
density by the replacement 

(3.7) 

e 

As a result, Eqs. (3.1) and (3.3) remain unchanged, while Eq.(3.2) is replaced by 

„,5/3 

y^HNJ) = W- 2 ^(Nt )a S /3 (ft ) 2 / 3 =pe 5 / 3 y(° ld \N,f,a ), (3.8) 

where a = TG^i, i.e. ao denotes what was denoted a st above. To keep track of the dependence on e (at a fixed G/i), 
we shall also rewrite Eq.(3.1) in terms of e and a , 

M/) ~ T-h^a^iftor^giy^Ml - m [a, /, W ncw) )]). (3-9) 

The function g[y] has the successive power-law behaviours g[y] oc y n with n = —1/3 for y < 1, n = —8/11 for 
1 < y < y eq , and n = —5/11 for y > y eq . In view of the scaling y( new ) oc pe 5//3 , Eq. (3.8), we see from Eq. (3.9) that 
hfi(f) will have successive power-law scalings with p and e of the form h^(f) oc e 2 / 3 (pe 5 / 3 )™ = p« e (2+5n)/3 Therefore 
the value of h^(f) for p < 1, e < 1 can be obtained from the corresponding value, for the same values of /, G/i, N and 
c, and for p = e = 1, by multiplying it with the following factors: p~ 1 / 3 e 1 / 9 for y < 1, p^ 8 / 11 e^ 6 / 11 for 1 < y < y eq , 
and p- 5 / 11 e~ 1 / 11 for y > y eq . 

The qualitative effect of varying p and e is now easy to understand. The graph of vs. G\x has a zigzag shape, 
with a rising line on the left, a short declining segment in the middle, and another rising line on the right. When 
p is lowered, at fixed Gfi and e, all three lines move up, with the central segment moving somewhat more than the 
right line, and the right line somewhat more than the left line. When e is lowered at constant G/i and p, the left line 
moves down, while the central segment and the right line move up. The displacements of the three lines in this case 
are very unequal. The displacement of the left and right lines is small, unless e changes by 5 orders of magnitude or 
more, while the displacement of the middle segment is quite noticeable, even if e is changed by only one or two orders 
of magnitude. It can be easily seen that the local maximum and minimum of the curve h^(Gfi) are both shifted to 
the right if e is decreased (see below) . When e and p are changed independently of G/j,, the slopes of all three lines 
remain unchanged. This would not be the case if, for example, e were a function of Gfi, as in Eq.(2.13). 

A complementary way to qualitatively grasp the effect of p < 1 and e < 1 on the plot of h^{f) (considered for a 
fixed value of /) versus G/j, is the following. In Eq. (3.9) above the dependence of h^(f) on G(i comes both through 

the prefactor ajy 3 and through the dependence of ?/ new ) upon ao appearing in Eq. (3.8), namely 
can, however, have a fast grasp at the location and height of the extrema of the zig-zagged graph of vs. G/z by 
using the inverse dependence ao oc (y( ncw )) 3 / 8 to reexpress as a function of y( ncw \ instead of ao, or Gfi = ao/50. 
Keeping all the factors, this leads to 

M/) - e- 3/ V 5/8 ^i(/, to, c, N)G(y^), (3.10) 
where we have introduced the auxiliary functions 

^10 2 c 



H 1 (f,t ,c,N) = ± 



Nt 



5/8 

(ft r 3/4 (3.11) 



and 
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FIG. 1. Effect of a reconnection probability 10 < p < 1 on the gravitational wave amplitude of bursts emitted by cosmic 
string cusps in the LIGO/VIRGO frequency band (fu go = 150 Hz), as a function of the string tension parameter G p (in a 
base-10 log-log plot). Here, as in the following figures, the average number of cusps per loop oscillation is assumed to be c = 1. 
The horizontal dashed lines indicate the one sigma noise levels (after optimal hltering) of LIGO 1 (initial detector) and LIGO 2 
(advanced configuration). 



G[y] ee y*/*g[y] = y 7 / 24 (l + y)- 13 / 33 (l + y/y eq f"\ (3.12) 

The function G(y) has a maximum at y ~ 1, and a minimum at y ~ y eq . These being numerically fixed values, the 
result Eq. (3.10) above immediately shows that the heights of the extrema of the plot vs. Gp, depend on e and p 
only through the simple prefactor in (3.10). Namely, the heights of the extrema are proportional to 

^extrema ^ £ -3/8 p -5/8^ (3^3) 

Interestingly, the effect of e < 1 or p < 1 is always to increase the extremal values of the plot h^ vs. Gp. As for the 
locations of these extrema on the Gp axis, they are given by inverting the relation (3.8) linking y to Gp. This yields 
(modulo factors depending only on c, N, f and to) 

Gp, cx e- 5 /V 3/8 y* /8 (3.14) 

where y» denotes the location of an extremum on the y axis ,i.e. y* ~ 1 or y» ~ y eq . These being numerically 
fixed values, we see that the locations on the Gp axis of the extrema of h^(f) vary with e and p simply through the 
prefactor e~ 5 / 8 p -3 / 8 in Eq. (3.14). This factor is always larger than 1 when e < 1 or p < 1. In summary, the effects 
of e < 1 and p < 1 on the zig-zagged graph of vs. Gn, is to move it up by a factor e~ 3 / 8 p~ 5//8 , and right by a 
factor e~ 5 / 8 p~ 3 / 8 . It is easily checked that this simple behaviour is compatible with the more detailed results above 
about the "motion" with e and p of the three lines making up the zig-zagged plot h^(Gfi), when taking into account 
the successive logarithmic slopes (+7/9, —3/11, +5/11 [17]) of h^(f) versus Gfi. 

How do these modifications affect the detectability of the GW bursts from cusps? A smaller value of p can only 
increase h^, and therefore improves the detectability. The improvement is moderate when considering a given value 
of Gfi: for p ~ 10~ 3 , we gain one order of magnitude in h. The improvements in the detectability of GW burst signals 
in the LIGO and LISA detectors brought by reducing the value of p (from 1 to 10~ 3 , in steps of 10 _1 ) is illustrated 
in Figs. 1 and 2. 

The horizontal dashed lines in these figures indicate the one sigma noise levels, after optimal filtering, for detecting 
GW bursts in the corresponding detectors. In the "LIGO" figures, the upper horizontal line corresponds to the initial 
sensitivity, LIGO 1, or equivalently VIRGO, while the lower line corresponds to the planned advanced configuration 
LIGO 2. For details of the signal to noise analysis, see [17]. 
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FIG. 2. Effect of a reconnection probability 10~ 3 < p < 1 on the gravitational wave amplitude of bursts emitted by cosmic 
string cusps in the LISA frequency band (fu sa = 3.88 x 10 Hz) , as a function of the string tension parameter G fi (in 
a base-10 log-log plot). The horizontal dashed line indicates the one sigma noise level (after optimal filtering) of the LISA 
detector. 



Turning to the effect of e < 1, we find that the burst amplitudes h^(f) are very weakly affected by a decrease of e 
by a few orders of magnitude. This follows from the fact that the limiting sensitivities of LIGO and LISA correspond 
to the regime y < 1, where the amplitude h is only lowered by a factor e 1 / 9 . From this factor it would seem that 
it is only in the case where e < 10~ 10 that detectability by LIGO or LISA will be significantly affected. However, 
when e gets that small the O-function factor in Eq. (3.9) starts playing an important role even at LIGO or LISA 
frequencies, especially when considering the types of values of the string tension, Gli ~ 10~ 10 , suggested by recent 
stringy implementations of brane inflation [26,24]. Indeed, the crucial numerical factor in the O-function cut-off is 
the product aft = eTGp,ft ~ 10 9 - 2 e(Gn/10- 10 )(f /Hz) entering 6 m , Eq. (3.5). Therefore, when G\i ~ 1CT 10 , the 
cut-off brought by the O function starts significantly affecting the LIGO signal (Jligo ~ 100 Hz) when e < 10 -11 , 
and the LISA one (Jlisa ~ 10~ 2 Hz) when e < 1(T 7 . The effect of e < 1 on the detectability of GW bursts by LIGO 
and LISA is illustrated in Figs. 3 and 4. 

Let us now briefly discuss the case e > 1. In that case, the effective value of a when a is used to parametrize the 
typical size of a loop at cosmic time t is a e ff = TG/i, corresponding to e e // = 1, while the loop density is generally 
given by Eq.(2.17). The latter result takes different explicit forms, according to whether one considers the cases of 
loops formed in the matter era (then one gets Eq.(2.18)), loops formed in the radiation era and decaying in matter 
era, or loops both formed and decaying in the radiation era (see [1]). For loops formed in the matter era, which 
correspond to the most relevant part of the graph in terms of observability, the loop density Eq.(2.18) coincides with 
the e = 1 limit of the case e < l,p < 1 considered above. Therefore, the plot of h vs. G[i is obtained from the 
discussion above by keeping the effect of p < 1, while setting e = 1. In that case, we see that the net effect is to 
increase the detectability of GW bursts. 

The overall conclusion is that the results of our previous work [16,17] are quite robust 2 against the inclusion of the 
modifications parametrized by p < 1 and by e (both when e < 1 and e > 1), at least when e > 1CT 11 in the case of 
LIGO, or e > 10~ 7 in the case of LISA. In particular, it is notable that a smaller reconnection probability can only 
increase the detectability of cosmic superstrings, and that even a very small e leads to a vast range of detectability 
for the GW bursts emitted by a network of strings. However, smaller values of e might lead to cutting off the burst 
signal in GW detectors, when the string tension is Gli ~ 10~ 10 or less. See Figs. 1-4. 



2 Let us also mention that Ref. [40] has shown that the waveform of GW bursts from cusps derived in [16,17] is robust against 
the presence of small-scale wiggles on the strings. 
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FIG. 3. Effect of a smaller fractional loop-length parameter 1CP 12 < e = a/(50G/i) < 1 on the gravitational wave amplitude 
of bursts emitted by cosmic string cusps in the LIGO/VIRGO frequency band (fu go = 150 Hz), as a function of the string 
tension parameter G \x (in a base-10 log-log plot). The horizontal dashed lines indicate the one sigma noise levels (after optimal 
filtering) of LIGO 1 (initial detector) and LIGO 2 (advanced configuration). 




FIG. 4. Effect of a smaller fractional loop-length parameter 10~ < e = a/(50G/i) < 1 on the gravitational wave amplitude 
of bursts emitted by cosmic string cusps in the LISA frequency band (fu sa = 3.88 x 1CF 3 Hz) , as a function of the string 
tension parameter G fi (in a base-10 log-log plot). The horizontal dashed line indicates the one sigma noise level (after optimal 
filtering) of the LISA detector. 
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IV. STOCHASTIC GRAVITATIONAL WAVE BACKGROUND 



We now consider the stochastic GW background produced by oscillating string loops. It was shown in [16,17] 
that this background appears, in general, as the superposition of occasional (non Gaussian) bursts, on top of a 
nearly Gaussian "confusion noise" h 2 on f usion (f), made of overlapping bursts. For some (rather large) values of the 
effective loop-length parameter a e //, Eq. (2.19), namely a e ff <~ 10~ 5 in the standard case, the individual burst 
events detectable during a typical pulsar-timing observation time scale T ~ 10 yr have an amplitude which might be 
comparable to the confusion background. This raises subtle issues about the detectability of such a mixed Gaussian- 
non-Gaussian background. In the discussion of this section, we shall, however, restrict ourselves to considering the 
simpler case of relatively small values of the effective loop-length parameter where the individual, non-overlapping, 
bursts are negligible compared to the background "confusion noise". Note, also, that, even in this case, Refs. [16,17] 
found that the quantity usually considered in the cosmic string literature, namely the "rms" noise h 2 ms (f), averaged 
over all bursts, was not necessarily a good estimate of the observationally relevant confusion noise, h 2 on j usion {f). 
Indeed, h 2 ms (f) includes, contrary to h 2 on j.- usion (f), the time-average contribution of rare, intense bursts, which are 
not relevant to a pulsar experiment of limited duration. 

The confusion noise can be written as the following integral over the redshift z (using Eq. (6.17) of [17] with the 
replacement Eq. (3.6) above) 

/r/z 
-n(f,z)h 2 (f,z)e[n(f,z)-l], (4.1) 

where 

n(f, z) = 10 2 ^(/t )- 5 / 3 a - 8 /V^)C(z), (4.2) 

is the number of cusp events (in a time window ~ / _1 and around redshift z) generating GW bursts around frequency 
/. Here 

M*) = + ^r 7/6 (! + ^e 9 ) 11/6 , (4.3) 



C[z] = l + 9z/(z + z eq ), (4.4) 

and 

h(f,z) = G m 2 / 3 (ft )- 1 / 3 Vh (z)e(l-e m [a,f,z}), (4.5) 

with 

<p h (z) = z- l (l + z)- 1 /3(i + z/z^)- 1 / 3 . (4.6) 

The factor C[z] interpolates between 1 in the matter era and 10 in the radiation era. It is incorporated to refine, in 
the standard case at least, the order of magnitude estimate Eq. (2.14) of the loop density. The quantity 9 m [a, /, z] 
entering the step function in the burst amplitude h(f, z) is that defined in Eq.(3.5) above. As above the step function 
0(1 — m [a, f, z]) cuts off the Fourier components that would correspond to a mode number m < 1. Note the presence 
of a second step function, Q[n(f,z) — 1] in h\ on f usion {f), which (approximately) limits the integral over the redshift 
z to the overlapping bursts. Indeed, the "confusion noise" Eq.(4.1) is obtained by subtracting from the "rms noise" 

h 2 rms (f) = Jv(dz/z)n(f,z)h 2 (f,z) the signals of low redshift jf f] \dz / ' z)n(f \z)h 2 {f ', z) with n(f,z c (f)) ~ 1. As 
shown in [17], the latter signals can be intense, but they are rare and non-overlapping, their total number (in a time 

window ~ being ~ f* c ^\dz/z)n(f, z) ~ n(f, z c (f)) ~ 1. (We use here the fact that n(f, z) is a monotonically 
increasing, power-law-type, function of z.) 

Finally, we associate to the dimensionless squared GW amplitude (per octave of frequency) h 2 con f usion (f) its energy 
density per octave of frequency, fdp g /df - [2tt fh confusion {f )} 2 /(lQirG) ~ (w / ' ^G)f 2 h 2 confusion {f), and the corre- 
sponding fractional contribution (per octave of frequency) of the confusion GW noise to the cosmological closure 
density, Q g (f) = (f / p c )(dp g /df), where p c - l/(67rGig) (see Eq. (6.20) of [17]) 

3-7T 2 

n g ° n f—(f) ~ —(ft ) 2 h 2 confusion (f)- (4-7) 
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FIG. 5. Effect of a reconnect ion probability 10 -3 < p < 1 on the fractional contribution fl g (f pBr ) (around the frequency 
fpsr ~ l/(10yr)J to the cosmological closure density of the stochastic GW noise due to overlapping GW bursts emitted by 
a network of strings. [Base-W log-log plot.] The upper (solid) horizontal line indicates the upper limit Q. g < 6h~ 2 x 10~ 8 
derived from 8 years of high-precision timing of two millisecond pulsars: PSR 1855+09 and PSR 1937+21. The middle (dashed) 
horizontal line indicates the potential sensitivity of 17 years of high-precision timing of PSR 1855+09 (see text). The lower 
(dashed) horizontal line indicates the expected sensitivity from the timing of the set of pulsars to be hopefully detected by a 
square-kilometer-array of radio telescopes. 



This quantity is plotted, for a pulsar timing frequency f psr ~ O.lyr -1 corresponding to a typical ~ 10 yr obser- 
vational window, as a function of Gfi, and for various values of p and e, in Figs. 5 and 6. The horizontal lines 
in these figures correspond to various realized, or planned, pulsar timing experiments. The upper (continuous) line 
corresponds to the (95% confidence level) upper limit Q g h 2 < 6 x 1CP 8 derived in [14] from 8 years of high-precision 
timing of two millisecond pulsars: PSR 1855+09 and PSR 1937+21. [Note that a Bayesian reanalysis of the data of 
[14] gave, under the choice of "Jeffrey's prior", the slightly less stringent limit Vt g h 2 < 9.3 x 10~ 8 [15]. Note also that, 
consistently with our use of H ~ 65 km/s/Mpc and t ~ 1.0 x 10 10 yr, we have h 2 ~ (65/100) 2 ~ 0.42] 

Recently, the data set for these two pulsars has been extended to a 17- year continuous span by piecing together 
data obtained from three different observing projects [34]. This is the first realization of the concept of Pulsar Timing 
Array (PTA) [41,42]. However, the upper limit on £l g h 2 that one can deduce from this 17-year combined data set is 
unclear to us. On the one hand, [34] computes two widely different upper limits by using two different approaches. 
A Neyman-Pearson test leads, according to [34], to a 95% confidence level limit of only £l g h 2 < 2.8 x 10 -6 , which is 
much less stringent than the limit obtained in [14] from 8 years of data. In view of this surprising result, Ref. [34] 
then resorted to a rather coarse estimate of an upper limit, based only on saying that "the largest amplitude sinusoid 
that one could conceivably fit to the PSR 1855+09 data" is 3 fis, for a frequency / = 1/(17) yr. This led to the 
upper limit Q, g h 2 < 2 x 10~ 9 , which is now more than ten times more stringent than the limit based on 8 years of 
data. On the other hand, a look at the PSR 1855+09 residuals reported in [34] shows that, during an intermediate 
period of ~ 5yr corresponding to the Green Bank data, there were residuals reaching the ~ 30/zs level, i.e. a much 
larger level than the ~ 3/zs level typical of the (pre- and post-upgrade) Arecibo data. This large "activity" of the 
pulsar data points during the Green-Bank-only period is responsible for the Q, g h 2 < 2.8 x 10~ 6 Neyman-Pearson-test 
limit. Apparently, the author of [34] seems to favor the tighter limit Q, g h 2 < 2 x 10~ 9 , obtained without using any 
statistical reasoning, more than the Neyman-Pearson-test one Vl g h 2 < 2.8 x 10~ 6 . Personally, in view of the "large 
activity" exhibited by the Green Bank data, we are not convinced that the level Q g h 2 = 2 x 10 -9 can be considered 
as a real upper limit on £l g h 2 . Pending a more complete statistical analysis of the present 17-year combined pulsar 
data, we can only consider this level as the potential sensitivity level of such an extended data set. Accordingly, we 
have represented the level Q g h 2 = 2 x 10~ 9 as a dashed line in Figs. 5, 6. 

Actually, wc wish to suggest, with due reserve, that the large scatter recorded in the Green Bank data set might 
be due to the real effect of a transient GW burst activity of some sort, e.g. to a near (z <C 1) or rare (c -c 1) cusp 
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FIG. 6. Effect of a smaller fractional loop-length parameter 10 -3 < e = a/(50Gfx) < 1 on the fractional contribution Q g (f psr ) 
(around the frequency f par ~ l/(10yr) ) to the cosmological closure density of the stochastic GW noise due to overlapping GW 
bursts emitted by a network of strings. [Base-10 log-log plot.] The upper (solid) horizontal line indicates the upper limit 
fig < 6h~ 2 x 1(T 8 derived from 8 years of high-precision timing of two millisecond pulsars: PSR 1855+09 and PSR 1937+21. 
The middle (dashed) horizontal line indicates the potential sensitivity of 17 years of high-precision timing of PSR 1855+09 (see 
text). The lower (dashed) horizontal line indicates the expected sensitivity from the timing of the set of pulsars to be hopefully 
detected by a square-kilometer-array of radio telescopes. 



event (indeed, there is no evidence of a steady stochastic red noise in the Arecibo data). When looking not only at 
the PSR 1855+09 data, but also at the PSR 1937+21 ones, one notices that, after fitting out a cubic term a P in 
the residuals, there remains some larger-than-usual activity, at the 4/xs level. It is clearly too early to use such data 
to fit for possible string-loop parameters, but we urge the observers not to dismiss this interesting possibility. 

Finally, looking ahead to the realization of the project of the Square Kilometer Array radio telescope [43], and 
of its consequent pulsar timing array, one can ultimately hope to reach, through pulsar timing, the sensitivity level 
Sigh 2 = 1CT 12 - 6 at a frequency f psr ~ O.lyr -1 [43]. This ultimate sensitivity level is indicated as the lowest dashed 
line in Figs. 5, 6. 

Figs. 5, 6 illustrate the effects of p and e on the detectability of the stochastic GW background generated by a 
cosmological network of (supcr)strings. This background is essentially made of the superposition of overlapping burst 
signals. Therefore we expect that the effects of p and e discussed above on the case of individual bursts will somehow 
extend to this "confusion noise" . Indeed, Fig. 5 shows that decreasing p increase the signal, and therefore improves 
its detectability by pulsar timing experiments. Concerning the effect of decreasing e the situation is a bit different 
from what happened in the case of individual bursts in the LIGO or LISA frequency band. The good news is that 
the amplitude of the stochastic background, versus Gli, in the extended domain where it is a rather flat rising line, 
increases when e decreases. The bad news is that, because of the low value of the pulsar frequency band f psr <~ O.lyr -1 , 
the stochastic signal is quite sensitive to the left cut-off brought by the step function 0(1 — 9 m [a, /, z]) present in each 
burst signal Eq. (4.5). This step function corresponds to saying that the Fourier series representing the GW amplitude 
emitted by a periodically oscillating loop only contains mode numbers |m| ~ (0 m (f, z))~ 3 > 1. As in the case of 
individual bursts considered above, the main numerical factor (when z <~ 1, which is the case on the left of the graph 
plotting fig vs. G^jl) which determines this cut-off is the product afto. When this product gets smaller than 1, the 
6 function cuts off the signal. Numerically, one has, for the pulsar frequency band, af psr to ~ 10 9 a <~ 5e(G^/10~ 10 ). 
Therefore, when considering, for instance, the type of values G/i ~ 10 -10 expected from brane inflation models 
[20,21,26,24], we see that values of e < 10 _1 are sufficient for cutting off the signal in the pulsar-timing band. On the 
other hand, for larger values of the string tension, the stochastic signal will instead increase when e decreases. 

In order to control analytically the values of the string tension that could be detected in pulsar timing exper- 
iments, it is useful to derive an approximate analytical approximation for the stochastic signal £i c ° n f uslon ( f:^ ^ 
(37r 2 /2)(/i ) 2 ^confMsiori(/)- By looking at the integrand of Eq. (4.1), one can see that in the low-frequency part of the 
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spectrum which is relevant for the millisecond pulsar observations (gently rising lines in Figs. 5, 6), this background 
is produced by the loops radiating at redshift z ~ 1, i.e. within the present Hubble radius, t ~ to- The value of the 
integral Eq. (4.1) can then be approximately estimated 3 by replacing the integral / dz/z by the value of the integrand 
at z ~ 1. This leads to 

1 n 2 

h 2 confusion (f) ~ -^cp-^-^G^YG^{ft Q )-V\ (4.8) 



and therefore, using Eq. (4.7), to 



Sl g {f) ~ SOG^cp-h-'/^TG^fto)- 1 / 3 , (4.9) 



where the numerical factor 30 comes from the combination 37r 2 10 2 /(2r). Note that, as was exhibited in Figs. 5, 6, a 
decrease in either pore from their standard model values (p = e = 1) increases the intensity of the background. Let 
us recall that these results are only valid on the gently rising slope on the left of Figs. 5,6, and in particular that 
they cannot be applied in the domain aftg < 1, corresponding to the left cut-off apparent in the figures. As said 
above, at frequencies relevant for millisecond pulsar measurements, / ~ 0.1 yr _1 , this cut off sets in when a < 10 -9 . 

For completeness, let us also sketch a direct derivation of the result (4.9) based on keeping track of the total energy 
emitted by the network of strings. A loop of length I radiates at a discrete set of frequencies, f m = 2m/ 1, but for 
large enough m it is well approximated by the continuous spectrum 

dP g /df ~ TG^l{.fl)- A '\ (4.10) 

The slow fall-off with frequency oc / ~ 4 / 3 holds when a cusp forms during a loop oscillation. For order of magnitude 
estimates, this formula can be used even for m ~ 1, but one has to remember that the spectrum is cut off at / < I /I. 

Assuming first that a <C TG/i, the total number of loops produced in our Hubble volume per Hubble time is 
N ~ I /{pa). Each loop radiates for a time period t <~ (a/TGp)t , so the GW energy density from all loops (with 
cusps) that radiated during the present Hubble time is (per octave of frequency) 



dp g dP g c pc 1/3 

f W~ f ^ T P<~p4 ifato) ■ (4 - n) 



Here, the factor c is added because it measures the fraction of the loops exhibiting cusps. Expressing this in units of 
the critical density, p c ~ l/(67rG^o), we get for fl g (f) = (/ / p c )dp g /df the result 

n fl (/) ~ GnGpcp-h-^iTGpfto)- 1 / 3 , (4.12) 

which is in good agreement (modulo a factor ~ 0.63) with the estimate (4.9) above. Let us recall that all our estimates 
try to keep the possibly important powers of 27r, but neglect various factors "of order 2" . [We note in this respect 
that r ~ 50 comes essentially from a factor (27r) 2 , while the factor 10 2 in Eq. (4.2) came from a factor oc 547T.] 
For completeness, let us mention that the case of a 3> TGp has been reviewed in [1]; the result is 

Sl g (f) ~ GirGpcp-^faoto)- 1 / 3 , (4.13) 

where «o = FGp, as above. 

Consistently with what we said above, the result (4.13), corresponding to e > 1, can be obtained from the result 
(4.9), corresponding to e < 1, by replacing e — ► e e // = 1 in the latter result. In order to treat both cases together we 
can therefore replace everywhere e by, say, e e // = e/(l + e). With this notation, our approximate analytical result 
(4.9) numerically yields 

n g (f)h* ~ W-^c(Gp)^p-h; f f(f/f psr )-^. (4.14) 
Alternatively, any pulsar timing sensitivity level fl g (f)h 2 corresponds to a bound on the string tension of order 

Gp ~ 10 3 - 7 (^( / )/, 2 ) 3 / 2 ( / // psr ) 1 / 2 c -3/ 2p 3/ 2e l/2^ (4 15) 



3 Note, however, that this estimate neglects a subdominant, but significant "floor" contribution coming from the radiation era, 
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Let us recall that c < 1 denotes the number of cusps occurring per loop oscillation. The parameter c is expected to 
be ~ 1 for generic smooth loops [33] , but the presence of many "kink" discontinuities along the loop might decrease 
the (effective) value of c below 1. It seems, however, reasonable to assume that c > 0.1. 

From Eq. (4.15), taking into account that in all cases we have p 3 ^ 2 e 1 Jj.j < 1, we get the inequality c 3 / 2 G/i < 
10 3 - 7 (O s (/)/i 2 ) 3 / 2 (/// psr ) 1 / 2 . Therefore the firm upper limit Q g h 2 < 6 x 1CT 8 [14] obtained for a frequency / ~ 
l/(7yr), yields the upper limit (?l 2 G[i < 10~ 7 , already quoted in the Introduction. Let us also consider the potential 
detectability ranges of pulsar timing arrays. The sensitivity level of the current realization of the pulsar timing array 
(PTA), n g h 2 = 2 x 1(T 9 for / - 1/(15. 6yr) (see above), corresponds to G^i ~ 3.6 x 10- 10 c- 3 / 2 p 3 / 2 e 1 J f 2 f . This is 
an impressive number which shows the vast discovery potential of pulsar timing experiments, and, in particular, the 
possibility for PTA experiments to detect the type of string tensions expected from recent superstring cosmology 
models [26,24]. One must, however, keep in mind the caveat exhibited in Fig. 6 about the adverse cut-off effect of 
having e < 1 to which pulsar experiments are especially sensitive. Fig. 6 shows that more sensitive PTA experiments, 
such as the ultimate Square Kilometer Array PTA, able to probe Q, g h 2 = 1CP 12 ' 6 , will become limited to the level 
G\i ~2x 10 _11 e _1 , by the ability of pulsar timing to probe the frequencies emitted by string loops. 

V. CONCLUSIONS 

It has been argued [20-22,26,24] that F- and D-string networks can naturally be formed at the end of brane inflation 
[23], with string tensions in the range 10~ n < G[i < 10~ 6 . We focussed on models where only one type of string 
(F or D) is formed, and considered the gravitational wave (GW) signatures of cosmological networks of such strings. 
We studied how the finding [16,17] that GW bursts emitted from cusps of oscillating loops should be detectable by 
LIGO and LISA interferometers for values of G\i as small as 10~ 13 might be modified by two separate effects. First, 
the reconnection probability p for intersecting F or D strings might be, contrary to the case of ordinary field-theory 
strings, significantly smaller than 1 [22,21,25]. Second, Refs. [16,17] had assumed that the characteristic size of newly 
formed loops was l(t) <~ at, with a ~ 50G/Z, as expected from from standard GW radiation- reaction arguments [11]. 
However, recent analyses [31,32] have suggested that gravitational radiation is less efficient than originally thought, 
and might result in a much smaller typical size for newly formed loops: l(t) ~ at, where a ~ e50G^, with f«l. 

A detailed analysis of the effects of the two parameters p and e on on the detectability of GW bursts by LIGO or 
LISA has shown that the results of [16,17] are quite robust, at least when e > 10 -11 in the case of LIGO, or e > 10~ 7 
in the case of LISA. See Figs. 1-4. In particular, it is notable that a smaller reconnection probability can only increase 
the detectability of cosmic superstrings. However, very small values of e might lead to cutting off the burst signals in 
GW detectors when the string tension is Gp ~ 10~ 10 or less. 

We have also considered the detectability, via pulsar timing observations, of the stochastic GW background produced 
by oscillating string loops. When the loop-length parameter a is large enough, namely a — 50eG/i > 10~ 9 , for the 
GW frequencies emitted by string loops at redshift z ~ 1 to fall within the pulsar sensitivity band / ~ l/(10yr), 
we find that the intensity of the GW background increases oc p _1 e -1 / 3 , when either p or e gets smaller than 1. In 
addition, we find that present pulsar timing experiments have the potential of detecting string tensions as small as 
Gfi <~ 3.6 x 10~ 10 c~ 3 ' 2 p s ' 2 e 1 ' 2 , where c denotes the number of cusps per string oscillation. We urge pulsar observers to 
reanalyze a recently obtained 17- year combined data set to see whether the large scatter exhibited by a fraction of the 
data might be due to a transient GW burst activity of some sort. We note that future versions of the "pulsar timing 
array" [34,41-43] might further improve the sensitivity of pulsar observations. The ultimate sensitivity of pulsar 
timing experiments might then become limited (by the ability of pulsar timing to probe the frequencies emitted by 
string loops) to the level G/j, ~ 2 x 10 _11 e _1 . In other words, if the suggestion of [31,32] is confirmed, and the reduced 
loop- length parameter e = a/(50G/i) turns out to take very small values e<<l (possibly of the form e ~ (50G/i) m , 
with m > 0), and if G/j, itself takes small values, G/i ~ 10~ 10 [26,24], the GW frequency band probed by pulsar timing 
experiments might fall in the domain a = 50eG/i < 10~ 9 where the GW spectrum from loops is cut off. In such a 
case, only higher- frequency GW experiments, such as LISA or, even better, LIGO, might be able to detect GW from 
string loops. 

Our conclusions show that it is urgent to develop a new generation of string network simulations able to determine 
how the crucial loop-length parameter a depends on GyU and p. It would also be quite important to determine the 
average number c of cusp events per loop oscillation. It is only when a and c are known that it will be possible to 
discuss with any reliability the detectability of cosmic superstrings by GW experiments. 
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